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We study the geodesic equation in the space-time of a Kerr black hole pierced by an infinitely 
thin cosmic string and give the complete set of analytical solutions of this equation for massive 
and massless particles in terms of Mino time that allows to decouple the r- and ^-component of the 
geodesic equation. The solutions of the geodesic equation can be classified according to the particle's 
energy and angular momentum, the mass and angular momentum per mass of the black hole. We 
give examples of orbits showing the influence of the cosmic string. We also discuss the perihelion 
shift and the Lense-Thirring effect for bound orbits and show that the presence of a cosmic string 
enhances both effects. Comparing our results with experimental data from the LAGEOS satellites 
we find an upper bound on the energy per unit length of a string piercing the earth which is 
approximately 10^^ kg/m. Our work has also applications to the recently suggested explanation of 
the alignment of the polarization vector of quasars using remnants of cosmic string decay in the 
form of primordial magnetic field loops. 

PACS numbers: 04.20.Jb, 02.30.Hq 



I. INTRODUCTION 



The motion of test particles (both massive and massless) provides the only experimentally feasible way to study 
the gravitational fields of objects such as black holes. Predictions about observable effects (light-like deflection, 
gravitational time-delay, perihelion shift and Lense-Thirring effect) can be made and compared with observations. 
Geodesies in black hole space-times in 4 dimensional Schwarzschild space-time [l| and Kerr and Kerr-Newman space- 
time have been discussed extensively. This has been extended to the cases of Schwarzschild-de Sitter space-times 
as well as to spherically symmetric higher dimensional space-times The analytical solutions of the geodesic 
equation in the Kerr space-time have been presented in |2, 5]. Moreover, solutions of the geodesic equation in the 
Kerr space-time have been given using elliptic functions 6] , while solutions representing bound orbits parameterized 
in terms of Mino time Q have been presented in 8]. Spherical orbits in the Kerr-(Anti)-de Sitter space-time have 
been discussed in @, while the general solution to the geodesic equation in (4-dimensional) Kerr-de Sitter Q and 
even general Plebanski-Demianski space-times without acceleration has been found [lo| . 

In a Fourier expansion has been used to compute the fundamental frequencies for bound orbits using Mino 
time. These results have direct application in the computation of gravitational waves that are created in extreme 
mass ratio inspirals, i.e. in binaries in which a stellar object moves on a bound orbit around a supermassive black 
hole. 

Gosmic strings have gained a lot of renewed interest over the past years due to their possible connection to string 
theory [l^ . These are topological defects [13i] that could have formed in one of the numerous phase transitions in the 
early universe due to the Kibble mechanism. Inflationary models resulting from string theory (e.g. brane inflation) 
predict the formation of cosmic string networks at the end of inflation fl^ . 

Different space-times containing cosmic strings have been discussed in the past. This study has mainly been 
motivated by the pioneering work of Bach and Weyl [l^ describing a pair of black holes held apart by an infinitely 
thin strut. This solution has later been reinterpreted in terms of cosmic strings describing a pair of black holes held 
apart by two cosmic strings extending to infinity in opposite direction. Consequently, a cosmic string piercing a 
Schwarzschild black hole has also been discussed, both in the thin string limit jl6j - where an analytic solution can 
be given - as well as using the full U(l) Abelian-Higgs model [13, 13: where only numerical solutions are available. 
In the latter case, these solutions have been interpreted to represent black hole solutions with long range "hair" and 
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are thus counterexamples to the No hair conjecture which states that black holes are uniquely characterized by their 
mass, charge and angular momentum. Interestingly, the solution found in [T6| is a Schwarzschild solution which 
however differs from the standard spherically symmetric case by the replacement of the angular variable by 
where the parameter /? is related to the deficit angle hy 6 — 27r(l — (3). In this sense, the space-time is thus not 
uniquely determined by the mass, but is described by the mass and deficit angle parameter /?. Of course, it is then 
easy to extend known analytically given space-times to the conical case by making the mentioned substitution. 

In order to understand details of gravitational fields of massive objects and to be able to predict observational 
consequences, it is important to understand how test particles move in these space-times. Next to the above mentioned 
examples the complete set of solutions to the geodesic equation in the space-time of a Schwarzschild black hole pierced 
by a cosmic string has been given recently in . 

The geodesic equation for the motion of test particles in the space-time of an uncharged, rotating black hole in 4 
space-time dimensions (given by the Kerr solution) pierced by an infinitely thin string aligned with the rotation axis 
of the black hole has been given in and spherical orbits and the Lense-Thirring precession have been discussed. 
Solutions to the geodesic equation in this space-time have also been given in [HI, but the test particle motion has 
been restricted to the equatorial plane and gravitomagnetic effects have been studied. Moreover, small perturbations 
around circular orbits have been discussed in [22]. The motion of a test scalar quantum particle in the space-time of a 
Kerr-Newman black hole pierced by a cosmic string has been discussed in [23| and it was observed that the presence 
of a cosmic string alters the corresponding observables. 

The aim of this paper is to determine the complete set of analytic solutions of the geodesic equations in the space- 
time of a Kerr black hole pierced by a cosmic string and to derive analytical expressions for observable effects which 
can be used for astrophysical searches for such cosmic strings. The accurate computation of geodesies for massive 
particles is important in order to understand gravitational wave signals from binaries which can later be compared 
with eventual gravitational wave measurements. The computation of geodesies for massless particles is important 
in order to understand how light signals pass by black holes or other massive objects. Moreover, our work has a 
direct link to the recently proposed explanation [2J| of the observed alignment of the polarization vector of quasars on 
cosmological scales [2^ . In [23| the assumption is used that two originally linked electroweak string loops decayed via 
the formation of monopole-antimonopole pairs in the early universe. The remnants of this decay are interconnected 
loops of magnetic field whose radii have grown due to the expansion of the universe and today should be on the order 
of Gpc. Interestingly, it was found that the rotation axis of a quasar would align with the direction of the magnetic 
field. Since the size of the magnetic field loops are much larger than the size of the quasars, we can assume the loop 
to be approximated by a straight line of magnetic field - that is aligned with the rotation axis of the supermassive 
black hole in the center of the quasar - at the position of the quasar. If we assume our infinitely thin cosmic string 
to be a toy model for finite width strings (e.g. electroweak or Abelian-Higgs strings) where the latter would have a 
magnetic flux along their axis, the model studied in this paper would describe a quasar with its rotation axis equal 
to the axis of the magnetic flux. 

Our paper is organized as follows: in Section II, we give the geodesic equation, in Section III, we classify the 
solutions, while in Section IV, we give the solutions to the geodesic equation and in Section V we present examples of 
orbits. In Section VI we discuss observables such as the perihelion shift and the Lense-Thirring effect. We conclude 
in Section VII. 



II. THE GEODESIC EQUATION 



We consider the geodesic equation 



d^x^ „„ dxP dx'^ 



, = (1) 

dT^ ' dT dT ' 

where F^^ denotes the Christoffel symbol given by 

^^p. = Ig^" {dp9.. + - a.ffp.) (2) 

and r is an affine parameter such that for time-like geodesies dr^ = g^jjdx'^dx'^ corresponds to proper time. 
The explicit form of the metric that we are studying in this paper reads 

dt + ^dr^ + p^dO^ + 13^ [r^ + sin^ ed(j? - P dtdcj) , (3 

p^ J A \ J 
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where — r"^ + a? cos^ 9, A = — 2Mr + a^, a = J/M is the angular momentum J per mass M of the black hole 
and < /? < 1 is the deficit angle parameter that is related to the deficit angle <5 = 27r(l — /?). This metric describes 
a Kerr black hole pierced by an infinitely thin cosmic string that is aligned with the rotation axis of the black hole. 
The deficit angle appears due to the presence of the cosmic string and can be expressed in terms of the energy per 
unit length p of the cosmic string : S — SirGfi ~ 87r(77/Mpi)^, where i] is the typical symmetry breaking scale at which 
the cosmic string formed and Mpi is the Planck mass. Note that we are using units such that G = c = 1. 

For a = this metric reduces to the Schwarzschild solution pierced by an infinitely thin cosmic string fioj, while 
for /3 = 1 we recover the standard Kerr solution. 

Surfaces with A = correspond to horizons of the Kerr solution with horizon radius r± = M ± V — , where 
the + sign corresponds to the event horizon, while the — sign corresponds to the Cauchy horizon. Surfaces with 
2Mr = are the static limit and define a Killing horizon. The domain between the event horizon and the static limit 
is the ergosphere. In the following, we are only interested in non-extremal black hole solutions, i.e. in solutions with 
> a? . Note that the localization of the horizons and the static limit, respectively are not altered by the presence 
of the cosmic string. 

The Boyer-Lindquist coordinates (r , 9, 0) are related to cartesian coordinate {x, y, z)hy x — yfp^+a? sin 9 cos , 
y — \J ^ sin 9 sin and z = r cos 9. Hence r = corresponds to a disc with a deficit angle 5, while the physical 
singularity at r = 0, = 7r/2 is a ring with a deficit angle. It is then clear - in contrast to the a = limit - 
that negative values of r are allowed. When crossing r = to negative values of r one enters into another conical 
space-time that however possesses no horizons. The Penrose diagram of this space-time looks exactly like that of the 
standard Kerr space-time Q, however each point in the diagram would correspond to a sphere with deficit angle 5. 

The Lagrangian reads 

1 dx^" dx" 1 

where e — —\ for massive test particles and e = for massless test particles, respectively. The Killing vectors of this 
space-time are ^ and ^ . The constants of motion are the energy E and the angular momentum which are given 
by the generalized momenta pt and 

dC ■ ■ ( 2Mr\ ■ 2Mar 



Pt ^ -^^--tgu -(t)gt4,^ [I- ^^jt + p=^^ sir? 9 <t>=:E , (5) 

dL ■ ■ 2Mar , . , (r^ + a^)^ - Aa^ sin^ 6* , • 
P4> = —j^tgt^ + <l>g4>4> = -l3 —sm^9t + l3^- — sin^ 9 cj) =: (31^ . (6) 



The dot denotes the differentiation with respect to the affine parameter r. Therefore the presence of a cosmic string 
aligned with the rotation axis of the Kerr black hole decreases the magnitude of the generalized momenta due to 
the parameter (3 as compared to the standard Kerr space-time. In particular, for a given energy E a particle on the 
surface defined by gu — (the static limit) has larger (p as compared to the standard Kerr case. 

There is another constant of motion, namely the Carter constant K [26| that appears when separating the Hamilton- 
Jacobi equations 



95 _ 1 
a7 ~ 2 



9'"'id^S){d.S) , (7) 



which persists in the presence of a cosmic string. S denotes the principal function of Hamilton for which we make the 
following Ansatz 

S=^eT-Et + ^L,cl, + Srir) + Sg{9) , (8) 

where Sr{r) and Sg{9) are functions of r and 9 only, respectively. 

By differentiating S with respect to the constants e, E, and Q .= K — {L^ — aE)'^, where Q is the modified 
Carter constant, we find the geodesic equations 



p^r = ±/RW, (9) 

±Ve{0) , (10) 



pH 



aP(r) 

p^f3c^ = {L, csc^ 9^aE) + , (11) 

A(r) 

pH = a{L,^aEsmH)+^-!-^^^p^ (12) 
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with 

6(6*) = Q -cos"^ 9 {Ll csc^ e - {E"^ + e)) , (13) 

P{r) = E{r^ +a^) - L,a , (14) 

R{r) = P{rf -A{r){Q + {L,-aE)^ -er^) . (15) 

Introducing a new parameter, the so-called Mino time [7] given hy dX = the r- and 0-componcnt of the geodesic 



equation can be decoupled 



dr 

dX 
d9 
dX 



= ±VW) , (16) 

= ±VW). (17) 
1 f csc^ e-aEdO aPir) dr \ 



dX '^^ A{r)y^R(r)dXJ 

dt _ a{L^- aEsw? 9) de {r^ + a^)P{r) dr 
dX ~ ^ A{r)y/R{rj dX 



(18) 
(19) 



For a = these equations reduce to the equations of motion in a space-time of a Schwarzschild black hole pierced by 
a cosmic string [l3|, while for /3 = 1 we recover the geodesic equation in the Kerr space-time. As initial conditions 
we will choose the +- signs in ([T5)) and PT)) . 



III. CLASSIFICATION OF SOLUTIONS 



The classification of the solutions of the geodesic equations ([T6|) -([T9 | can be done with respect to the modified 
Carter constant Q, the mass of the black hole M, the angular momentum per unit mass of the black hole a as well 
as the energy E and angular momentum of the massive (e = —1) or massless (e — 0) test particle. Apparently, 
the deficit parameter /3 does not appear in the r, 6 and t component of the geodesic equation and will hence only 
influence the (/)-motion. 



A. r— motion 



In order to have solutions of the geodesic equation for r we have to require R{r) > 0. Therefore solutions of the 
r-component of the geodesic equation exist only for specific choices of E and Lz- R{r) will have either 4 real, 2 real 
and 2 complex or 4 complex zeros. On the boundaries between the domains in the E^-L^-plane corresponding to 
these three different possibilities R{r) necessarily has a double zero. In order to find these boundaries we make the 
Ansatz R{r) — {r — k)^ {{E"^ + e)r^ + pir + P2) where k, pi and p2 have to be determined. We find that for R{r) = 
we have the following parametric expressions for E and 

_ -2k^ -ieMn^ + {ea^)K~ KM 
^('*) ~ 7^= ) (20) 



eMK^ + {K-2a^-ea^)K^-{-iMK + 2,eMa^)K^ + {ea'^ + a^K)K + a'^MK 
^z(K) = j= , (21) 

where 



P^ = Ak^ + 8eMK^ + {AK - Aae)K^ - 8MKk^ + ia^Kn^ . (22) 

Note that r — k corresponds to a double zero of R{r) and hence represents spherical orbits. 

Examples for the polynomial R{r) are given in FigH] and Fig [21 In Fig|T] we show the different domains in the 
iJ^-Lj-plane for M = 1, a — 0.8 and K = A. The blue line corresponds to = 1, the green line to k s] — 00, 0] and 
the red line to k e [0,+oo[. The plots of R{r) in the domains Rl to R5 are also shown. For Rl R{r) possesses 4 
complex zeros, for R2 and R5 2 real and 2 complex zeros, while for R3 and R4 there are four real zeros. The following 
type of orbits are then possible ^] : 
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• Flyby orbit: the test particle starts from ±00, reaches a minimal value of |r| and flies back to r = ±cx3. There 
are two flyby orbits in R2 and R4, respectively. 

• Transit orbit: the test particle starts from ±00, crosses r = and continues to r = =Foo. There is a transit 
orbit in Rl. 

• Bound orbit: the test particle oscillates in an interval [r2,ri] where ri > r2- There are two bound orbits in 
R3 and one bound orbit in R4 and R5, respectively. 

• Spherical orbit: this is a special bound orbit with r constant, r = k is a spherical orbit that can be stable or 
unstable. 

In Figl5]we show how the domains Rl to R5 change when changing the Carter constant K (upper three flgures). 
Obviously R3 decreases in size when increasing K. We also give the change of Rl to R5 for changing angular 
momentum per unit mass a (lower three figures). In this case, the domain R4 is increasing for increasing a. 



B. 0— motion 



It is obvious from the form of the geodesic equations that we should have 8(6') > 0, i.e. 6'-motion is allowed only 
for those for which 0(9) > 0. This in turn means that when we fix Q and a only particular values of E and Lz are 
allowed. 

1. equatorial motion with 6 = tt/2 ioi massive and massless particles: in this case, it follows from that Q = 0, 



I.e. 



Lz^aE±^/K . (23) 



This is shown in Fig 3(a) and Fig |3(b) [ where the blue and the red line indicate = aE+\/K and = aE—^/K, 
respectively. 

2. polar motion with 9 — or 6 — n: obviously in this case, we have to choose — and Q > —a?{E'^ + e). 

3. motion for < 6* < tt with 6 ^ ti/2: in this case, we find that (see Appendix for details) 



for massive particles and 



^l {E-VW^e){aH + K) 

- 2 ea ^ ' 



for massless particles. This is shown in Fig 3(a) and Fig [3(b)[ where the green line indicates L 



\{E — \/E'^ — e) (a^e + K) / (ea). The allowed domain in the Lz-E^-p\ane is the one above this green line. 

We show an example of the different domains corresponding to zeros of &{0) in FigHjfor M = 1, a = 0.8 and K = 4. 
The function 9(6*) is plotted for the domains denoted Tl to T5. In domains Tl, T2 and T5 the polynomial &{0) 
has no real zeros, while it has four real zeros in domain T3 and two real zeros in domain T4, respectively. Clearly, 
only in the domains T3 and TA solutions to the geodesic equations exist. In T3 there are two domains in 9 for which 
Q{0) > 0. Here the particles oscillate between 0,„i,i and O^^x but cannot cross 9 = ^. In domain T4 the particle can 
cross the equatorial plane 9 = 7r/2. 



IV. ANALYTICAL SOLUTIONS OF THE GEODESIC EQUATION 



In the following, we will present the analytical solutions of the geodesic equation for the r, 9, cf) and t component. 
We will present these solutions using Mino time A. 
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A. 7 — motion 



The polynomial R{r) can be written in the form 

4 

R{r) = {E' + e)l[{r^r,) , (26) 

where the rj denote the zeros of R{r) with ri the largest real zero. Using the coordinate transformation r — 1/y + ri 
the polynomial R{r) can be transformed to a 3rd order polynomial in y using 

dX = = ^ (27) 

with 

bo = I, (28) 

bi = (-ra - r4 + 3ri - ra) , (29) 

b2 = (r2r4 - 2rir2 + rars + 3ri^ - 2rir3 - 2rir4 + r3r4) , (30) 

63 = (-r3ri^ + rir2r3 + rir2r4 + rir3r4 - r2ri^ + ri^ - r4ri^ - r2r3r4) . (31) 
Introducing the variable 



w.h .= |/5;(^ , o = -J| ,32, 



we can write the solution of ([TB)) as follows : 



r(A) = ^ +ri| . (33) 

(-(^ - ^0) + Cr^h^h j + a 



The integration constant is given by 



a = / , ^! . , (34) 



^4z3 - 32^ - h 



where zq denotes the value of z that corresponds to the initial radius and 



= ' = i^^? ■ ^''^ 



Depending on the sign of the discriminant D — §2 — 27.g| we have different values for Cr- In the following, we will 
denote the zeros of the 3rd order polynomial in the square zero of (|34p as ei > 62 > 63. 

1. D > : In this case R{r) has four real zeros which corresponds to domains i?3 and R4 in Fig[T] 

For the outer bound orbit in domain R3 with minimal radius r2 and maximal radius ri we choose zq = e\ 
such that Cr = K{IC) / \/e\ — §3 ~: a)i where K{IC) is the complete elliptic integral of the first kind and K, = 
\/e2 — — 63 is the modulus of the elliptic integral. For the inner bound orbit in domain R3 with minimal 

radius r4 and maximal radius rj, we choose zq = e.2 such that Cr — K{lC)/\/ei — e^+iK{lC') / \/ei — 63 =: 1I11+UJ2 
where /C' ~ \/l — ■ For the flyby orbit in R4 with initial minimal radius ri, we choose zq = 00 such that 
Cr — 0, while for the inner bound orbit in R4 with initial maximal radius r2 and minimal radius r3 we choose 
Zq — P'2 such that Cr = Wl + ti)2 ■ 

2. D < : In this case R{r) has two real zeros which corresponds to domains R2 and R5 in Fig[TJ In the following, 
we will denote the real zero of the 3rd order polynomial in the square zero of p4p as 62 , while the two complex 
zeros are denoted by ei and 63. 

For the flyby orbit in domain R2 with initial minimal radius ri we choose zq = 00 such that Cr = 0. For the 
bound orbit in domain R5 with initial maximal radius ri we also choose zq — 00, i.e. Cr = 0. 

3. D — : In this case R{r) has one real zero. The orbit in this case will be a bound spherical orbit. 

Note that there is also the possibility of four complex zeros. This however would simply correspond to a transit 
orbit with r = —00 — ^ +00. We will not discuss this case in detail in this paper. 
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B. 0— motion 

The solution of ([T7| as a function of Mino time A is given by 



eix) 



where 



arccos 



1-1/3 



± 



1 



yAip(^(-^-Ao)+Ce;g2,53) + 



(36) 



(37) 



The positive (negative) sign of the arccos corresponds to the choice 9 > 7r/2 {9 < 7r/2). For the 6 motion in domain 
T4, the solutions thus have to be "glued together" at 6' = 7r/2. 
The integration constant Cg is given by 



C'e 



dw 



92W - 53 



(38) 



where w is related to 9 by cos^ 9 = {^w + v) ^ , wq corresponds to the value of w that represents the initial 9 and 



.92 



22/3(a2)i/3(3a^a3_a: 



3ai 



93 



20,2 ~ 9aia2a3 
27af 



(39) 



where ai = 4Q, 02 = 4 (^a?{E'^ + e) — — Q) and 03 — —4a^{E'^ + e). The discriminant f = ~ 27g| is always 
positive, so we have three real zeros of the 3rd order polynomial, which we call ei > 62 > 63 in the following. Note that 
though we might have three real zeros for the polynomial in w these zeros might not fulfill {fiw + ly)^^ — cos^ 9 < 1. 
Moreover to each zero of the polynomial in w correspond two values of 9 fulfilling cos 61 — ±(/ii(; + v)^^^^. Hence, 
depending on the values of ei, 62 and 63 the polynomial Q{9) can have four, two or no real zeros. 

In T4, we typically choose wq = ei such that Cg = K{IC)/y/ei — 63 =: wi, where K{IC) is the complete elliptic 
integral of the first kind and K. = \/e2 — 63/ ^/ei — 63 is the modulus of the elliptic integral. 

In T3, we typically choose wq — 63 such that Cg — iK{JC')/^/ei — 63 —: W2, where K{IC') — \J\ — K?. 



-motion 



The equation for the component ([T8|) is separated into a 0-dependent and r-dependent part 

d9 H , ' dr =: dig + dir ■ 



The solutions for Ig and 1^ are (see Appendix for more details) : 

2 



= {L.-aE){X-X,)+Y."ArY 



(A - Ao) 



C,{xi) + In {a{x - Xi)) - In {a{xQ - xt)) 



Ir = -i^o(A-Ao)+ V - 

' ^ en' 



^1 P'("jO 



-(A- Ao) 



C(Mji) + In [a{u - Uji)) - In {a{uo - Uji)) 



(40) 

(41) 
(42) 



where C and a denote the Weierstrass zeta- and sigma-function, respectively, /i and v are the variables defined in p7p , 

do = {1 — v) / fj, and x = X/fi. Moreover, we have introduced the variable m = — ^i(A — Aq) + C,-^ with u{Xo) = —Cr 

and the 7 given in (P^ . In addition uji =ej with p(uii) = p(wi2) = ei and p(m2i) = p(w22) = 62- Finally, the Kj 
appear when rewriting dIr as follows 



dIr — Kq 



jdz 



dz 



y^Az^ - g2Z - 53 . ^ (z-ej)v/4z 



322: - 53 



(43) 



where e,- := — 



and 2/ = (?' — ri ) ^. 
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D. motion 

From ([T9)) we have 

rift dr - - 

dt = a{L, - aEsin^ 0) Z__ + (r^ + a^)A{r)-' P{r) " =: dig + dl,.. (44) 

WW) WW) 

To find analytical expressions for 1$ and Ir we proceed as in the case of Ig and Ir, respectively (see also Appendix) 
and find 



1 

Ig = a{L, - aE){\ - Aq) + a'sy^—^ 



— — —C,{xi) + In {a{x - Xi)) ~ In {a{xQ - x,)) 



(45) 



where p(xi) = p(a;2) = j [i. The simple poles of the function (j63p given by xi, a;2 are in the fundamental domain 
\2aijj\ + 26(u;2|a, 6 e [0, 1]} where 2a;i G R and 2a;2 G C. In addition, we have 



4 2 



/.^Co(A~Ao) + EE^ 



— ^C(^ji) + In (cr(w - Wji)) - In ((t(uo - Uji)) 



(46) 



where d are the coefficients from the partial fraction and and the poles of the rational function F{z) from the 
transformation (r^ + o^)A(r)~-^P(r) f""^ . In addition p('«ii) — p(u2i) = e^. 



V. EXAMPLES OF ORBITS 

In the following we give plots of geodesies of massive and massless particles, respectively in the space-time of a 
non-extremal Kerr black hole pierced by an infinitely thin cosmic string. In particular we will demonstrate how the 
presence of the cosmic string alters the test particle motion. 

A. Motion of massive test particles 

The domain of existence of solutions of the geodesic equation can be obtain from the intersection of the allowed 
domains of the L^-S^-plane obtained from the requirement &{0) > and R{r) > 0, respectively. This leads to four 
domains in the L^-iJ^-plane. These are shown for M = 1, a = 0.8 and K = 12 in Figl5]and denoted by Zl to Z4. 

1. Zl: This domain is the combination of domain TA from Fig|4]and domain _R3 from Fig[T] The possible orbits 
are two bound orbits on which the test particle can cross the equatorial plane at 6* = 7r/2. The effect of the 
cosmic string on the geodesic motion on the outer bound and inner bound orbit is shown in FigjS] and Fig|71 
respectively. Note that the test particle on the inner bound orbit (see FiglT]) crosses the event and the Cauchy 
horizon (red circles) several times. 

2. Z2: This domain is the combination of domain TA from Fig|4]and domain RA from Fig[T] The possible orbits 
are one flyby and one bound orbit on which the test particle can cross the equatorial plane at = 7r/2. The 
effect of the cosmic string on the flyby orbit is shown in FiglS) 

3. Z3: This domain is the combination of domain T4 from FigUand domain R2 from FigHJ 

The possible orbits are two flyby orbits (from ri to oo and — oo to r2) on which the test particle can cross the 
equatorial plane at = tt/2. The effect of the cosmic string on the flyby orbit (from ri to oo) is shown in FiglHl 
Note that the test particle crosses the event and the Cauchy horizon (red circles) several times. 

4. Z4: This domain is the combination of domain T4 from Fig|4]and domain R5 from Fig[TJ 

The possible orbit is one bound orbit on which the test particle can cross the equatorial plane a,t = 7r/2. The 
effect of the cosmic string on this bound orbit is shown in FigfTOl Note that the test particle crosses the event 
and the Cauchy horizon (red circles) several times. 
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B. Motion of massless test particles 



The domain of existence of solutions of the geodesic equation can be obtain from the intersection of the allowed 
domains of the L^-i?^-plane obtained from the requirement Q{9) > and R{r) > 0, respectively. This leads to three 
domains in the L^-iJ^-plane. These are shown for M = 1, a = 0.8 and X = 12 in Fig. [Til and denoted by Nl to N3. 

1. Nl: The possible orbit is a cross-over flyby orbit on which the test particle cannot cross the equatorial plane 
a.t 9 = tt/2. The effect of the cosmic string on the cross-over flyby orbit (from ri to oo) is shown in FigfTSl 

2. N2: The possible orbits are one bound and one flyby orbit on which the test particle can cross the equatorial 
plane at 9 = Tr/2. The effect of the cosmic string on the cross-over flyby orbit (from ri to oo) is shown in Fig|51 

3. N3: The possible orbits are one flyby and one bound orbit on which the test particle can cross the equatorial 
plane at 9 = 7r/2. The effect of the cosmic string on the flyby and on the bound orbit is shown in Fig 11^ and 
Figim respectively. 



VI. OBSERVABLES 



For the standard Kerr space-time {(3 — 1) the analytical expression of the so-called fundamental frequencies of 
bound orbits have been given in j7l| in the frequency domain using a Fourier transformation and in Q as a function 
of a Jacobi elliptic integral. Here, we will present the fundamental frequencies for /3 < 1 in the form of a Weierstrass 
elliptic integral p and analyze their dependence on the deficit angle. 

For bound orbits r(A) and 9{X) become periodic functions that are independent of each other. We can then define 
the fundamental frequencies of these bound orbits using Mino time A. Using the results of Sections [IV Al and IIV Bl we 
find 



A^ 



dr 



2"fLUl 



dcos9 
^e(cos6') 



(47) 



where r(A) = r(A + nAr) and 9{\) = 9{X + nkg) 
and the apoapsis in the radial direction, while 9-^ 
angular frequencies then read 



for any integer n [11|. rmin and rmax correspond to the periapsis 
in corresponds to the minimal value of the polar coordinate. The 



T - — - ^ 

Kr "fUJrl 



2tT TT 

= -r— = t; 

Ag 2flUJi 



(48) 



In Q the frequencies for the (p- and i-component have been defined via an average over the orbital periods A^ and 
Ag. Here we read these periods off from our solutions for the 0- and t-component of the geodesic equation. These 
solutions contain terms that correspond to oscillations with periods A^ and Ag , respectively and a term that describes 
a linear increase in Mino time A. and F are the frequencies of 4> and t in Mino time, respectively that correspond 
to this linear increase [1] . By using the results of sections IIV CI and IIVDI we can find the expressions for these 
frequencies. This gives 



= ^ \ {L,-aE-Ko)+ ^ 



(49) 



and 



F = aiL,-aE) + Co + Y,J2 

i=i ]=i 



a^E- 



1 C{Xi) d C{Uji) 



p'{xi) fi p'iuji) 7 



(50) 



These results differ from that of the standard Kerr space-time by the factor of 1//3 in 

As shown in ll| the angular frequencies calculated using Mino time A are related to the angular frequencies ilr, 
fig and fl^ calculated using a distant observer time as follows 



fig 



(51) 
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If these frequencies are different they give rise to the precession of the orbital elhpse and of the orbital plane. In 
particular, the perihelion shift is related to the difference between the angular frequency of the radial motion fi^ and 
the angular frequency of the (^-motion, while the Lense-Thirring (LT) precession is related to the difference between 
the frequencies of the two angular motions 

i^pcrihclion = — fir i ^LT = ^cf) — ■ (52) 

In comparison to the standard Kerr space-time with (3 = 1, the frequencies of the perihelion shift riperihciion and the 
Lense-Thirring precession J^lt are hence bigger for (3 < 1, i.e. in the Kerr space-time including a deficit angle. This 
is e.g. clearly seen in Fig. 9, where the perihelion shift increases for decreasing /3 (see plots in the x-y-plane) and the 
increase of precession is seen when studying the 3-d orbits. 

Using the results from the LAGEOS satellites [l^l we can estimate an upper bound for the deficit parameter and 
hence for the energy per unit length of a cosmic string piercing the earth. The theoretical value of the JIlt for the 
earth is given by 39 • 10~^ arcseconds/year. The LAGEOS satellites have measured this value with an accuracy of 
10%. Using this we find that 

^ltW ^l)- ^ltW = 1) = ( 3 - 1 ) ^0(/5 1) < 4 • lO-^arcseconds/year . (53) 



Assuming that Qtp{l3 = 1) is approximately 2tt per day, we find that 

4 - 1 < 10-" ^ 7^ < 10"" . (54) 

P ZTT 

This transfers to a bound on the energy per unit length /i of a cosmic string piercing the earth that reads 

fi < 10^6 kg/m . (55) 

Surprisingly, this upper bound is in very good agreement with that found for a cosmic string piercing a Schwarzschild 
black hole when comparing the theoretical results with the experimental results for the perihelion shift of Mercury 
and the light deflection by the Sun fl9| . 



VII. CONCLUSIONS 



It would be fascinating to detect cosmic strings in the universe since this detection would open the window to 
the physics of the very early universe and could prove (or disprove) theories such as string theory, supersymmetry 
and Grand Unification. In most cases the prediction of the detection of these objects has focused on the Cosmic 
Microwave background (CMB) data [2^. Here we discuss another possibility, namely that cosmic strings might be 
detected due to the way that test particles move in their space-time. Hence, in this paper we have studied the 
analytical solutions of the geodesic equation in the space-time of a Kerr black hole pierced by an infinitely thin cosmic 
string aligned with the rotation axis of the black hole. We have given the analytical solutions for the t-, r-, 9- and 
0-components of the geodesic equation in terms of Mino time which allows to decouple the r- and 0-motion. We see 
that the main difference to the standard Kerr space-time (which corresponds to the limit of vanishing deficit angle) is 
the change of the 0-motion. In particular the precession of the orbital plane as well as of the orbital ellipse of bound 
orbits will increase for increasing deficit angle, i.e. increasing energy per unit length of the string. Comparing our 
results with the LAGEOS data, we find that the upper limit for the energy per unit length of a cosmic string piercing 
the earth is fj, < 10^^ kg/m. Our results have also applications in the computation of gravitational wave templates for 
extreme mass ratio inspirals and to the recently suggested alignment of the polarization vector of quasars, respectively. 
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VIII. APPENDIX 



A. Domain of existence for the motion 



In order to find the expressions in p4|) and ((25)) we use the variable v :— cos^ 9 and rewrite ([T7| as follows 

dX = -^ (56) 

with Qy = — 4a^(i?^ +e)v^ +4 (^a^{E^ + e) ~ ~ Q) +4Qu. For v e]0 : 1[ we assume that two of the zeros merge, 
i.e. we can then write 0^ as follows 

Oy^iv- ^i)^^2V (57) 

where and ^2 can be given as functions of L^, E, e and Q. These relations can then be used to give Lz in terms of 
E. 



B. Computation of Ig 

Using the substitutions of section HVBI the expression dig reads : 

Lzdv aEdv ^{Lz — aE)dw Lz{v + ^do)dw 
die = 1 H 1=^ ~ — , : H : , (58) 

(l-l')v(0f) V0U V4w3-g2W-ff3 (w - do)v4w^ - - 33 

where do = ^77^- Let us define the new variable x and the constant xq by using the Weierstrass p function : 

and Xq :— f — j^^^^=^^^^= , (59) 



y/ 4w3 ~ - 53 Joe y/ Aw^ - g2W - 53 
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i.e. w = p(x; g2, gs). With this definition x can be solved as a function of Mino time A : 

A-Ao 



X — Xq 



In addition dw can be written as 



dw ~ p'{x)dx = ±\/4p'^(a;) — g2p{x) — g^dx where p 
Inserting dw into ((55)) gives : 



/ _ dp{x;g2,g3) 
dx 



die = ^i{Lz — aE)dx + {v + iido)- 



dx 



p{x) - do 

The first term in (j62p can be integrated easily. For the second term, we can use the formula 

2 



dx 



g p{x) -do ^ p'{xi) 



} s [{x - xo)({xi} + In {a{x - Xi}) - In {a{xo - Xi))] 



(60) 



(61) 



(62) 



(63) 



where p{xi) =p{x2) = do- The simple poles of the function in ()63p denoted by xi, X2 are in the fundamental domain 
{2auji + 2buj2\a, b £ [0, 1]}, where 2wi e K and 2a;2 G C. Hence the analytical solution of Ig is : 

d/e = le = fi(Lz - aE){x - xq) + ^ — — r- [{x - xo)C,{xi) + In ((7(2; - Xi)) - \n{a{xo - Xi))] 



2 



(A- Ao) 



C(a;j) + In {a{x - x^)) - In {(t{xo - x^)) 



(64) 



C. Computation of Ir 

With the substitution r = - + ri the expression dir can be written as a function of y as follows 



dir = 



aP{r)dr 



Moreover the integral in (j65p can be written in the form 



[Erl + Ea^ - L^a)y'^ + 2Eriy + E 



V 2ri-2M , ^ i 

y ^ r'f+a^~2Mri « ^ r^+a2-2Mri 



y-yi y-y2 



dy 
J Ru 



where 



2/1,2 



n- M± - a2 

-rf - a2 + 2Afri 



(65) 



(66) 



(67) 



and the Ki, i — 0,1,2 are constants. Using the transformation y = + a (where 7 and a are constants) the 
expression dir becomes : 



dir = Ko 



"fdz 



dz 



^4z3 - g2Z - gs {z - e^)^\z'^ - ~g2Z - .93 



where e, 



(68) 



We introduce the new variable u as follows 



13 



i.e. p{u; §2,93) = z. Then u can be given as a function of Mino time A. Finally using ([62)) - (|63)) the analytical solution 
of Ir is : 

K 

dir = Ir ^ iKq{u - Mo) + ^ ' . [{u " uo)C,{uji) + In (cr(u - Uj.i)) - In ((t(uo - 



-ifo(A-Ao)+ ^ - 



(A-Ao) 



Ciuji) + In (cr(u - Uji)) - In ((7(mo - Uji)) 



where m = — ^i(A — Ao) + C,.^ , uq = -C^ and Uji —Cj with p(mii) = p(ui2) = ei and p(u2i) = p(m22) = 62. 



(70) 




FIG. 1: The domain of existence of solutions to the geodesic equations for AI = 1, a = 0.8 and K = 4 that follows from the 
requirement R{r) > 0. The top left figure shows the domain of existence in the L^-E^-plane. The green and red lines represent 
R{r) — 0, where the green line is for k g] — oo : 0] and the red line for k £ [0 : +oo[. The remaining figures show the function 
R{r) for the different domains denoted by Rl to _R5. 



with a = 0.80 and K = 3.00 



with a = 0.80 and K = 8.00 



with a = 0.80 and K = 12.00 






(a) 



varying Carter constant K for AI = 1, a = 0.8. 



with a = 0.25 and K = 8.00 



with a = 0.60 and K = 8.00 




-20 
-25 
-30 

-35 - 



with a = 0.80 and K = 8.00 



(b) 



varying a for M = I, K = 3.0. 



FIG. 2: The change of the L^-E^ plot for varying Carter constant (a) and angular momentum per unit mass of the black hole 
a (b) respectively. 
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a = 0.80 and K = 3.00 



a = 0.80 and K = 6.00 



a = 0.80 and K = 12.00 




(a) 



M = 1, a = 0.8 with the change of K 



0.20 and K = 3.00 




a = 0.60 and K = 3.00 



a = 0.80 and K = 3.00 



20 40 60 80 100 120 





(b) 



a at M = 1, K = 3.0 with the change of a 



FIG. 3: The change of the domain of existence of solutions to the geodesic equations in the L^-E^ plane for several values of 
the Carter constant K and M — 1, a — 0.8 (a) and several values of the angular momentum per unit mass of the black hole a 



and M = 1 and K = 3|(b)[ respectively. Note that solutions exist only below the blue, above the red and above the green line. 



a = 0.80 and K = 4.00 Region 1 : Lz = -1 .78 and = 0.56 Region 2: Lz = -1 .38 and = 3.74 




0123 0123 0123 



FIG. 4: The domain of existence of solutions to the geodesic equations for M = 1, a = 0.8 and K = 4 that follows from the 
requirement 0(6') > 0. The top left figure shows the domain of existence in the L^^-E^ -plane . The remaining figures show the 
function Q{9) for the different domains denoted by Tl to T5. Solutions exist only for parameters Lz and E in domains T3 and 
T4, respectively. 
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with a = 0.8 and K = 12 




■4' ' ' ' ' ' ' ' 

0.5 1 1.5 2 2.5 3 3.5 



FIG. 5: The four domains Zl to Z4 in the Lz-E^-plane for which solutions to the geodesic equation for massive particles exist 
are shown for M — 1, a — 0.8 and K = 12. The red and blues lines come from the restriction Q{9) > and R{r) > 0, 
respectively. 




FIG. 6: The change of the outer bound orbit of a massive test particle (domain Zl) due to the change of the deficit parameter 
/?. Here Lz = 3.0, E = \/0.95, K — 12, M = 1 and a = 0.8. The red circles represent the radii of the event and Cauchy 
horizons, while the yellow circles denote the minimal and the maximal radius of the orbit, respectively. 
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FIG. 7: The change of the inner bound orbit of a massive test particle (domain Zl) due to the change of the deficit parameter 
/3. Here Lz — 3.0, E = \/0.95, K — 12, M = 1 and a — 0.8. The red circles represent the radii of the event and Cauchy 
horizons, while the yellow circles denote the minimal and the maximal radius of the orbit, respectively. 
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3D-flyby orbit 

r-e plane x-y plane 




(c) 13 = 0.50 

FIG. 8: The change of the flyby orbit (domain Z2) of a massive test particle due to the change of the deficit parameter /?. 
Here Lz — 3.104, E = 1.004, K = 12, M = 1 and a — 0.8. The red circles represent the radii of the event and Cauchy horizons, 
while the yellow circle denotes the minimal radius of the orbit. 




FIG. 9: The change of the flyby orbit (domain Z3) of a massive test particle due to the change of the deficit parameter /3. 
Here = —1, E — y/1.10, K = 12, M = 1 and a = 0.8. The red circles represent the radii of the event and Cauchy horizons, 
while the yellow circle denotes the minimal radius of the orbit. 
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3D-bound orbit 



r-9 plane 




x-y plane 




2 -2 



/3 = 1 



3D-bound orbit 



r-e plane 




x-y plane 
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3D-bound orbit 
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P = 0.70 



FIG. 10: The change of the bound orbit (domain Z4) of a massive test particle due to the change of the deficit parameter /3. 
Here — —1, E — \/ 0.5, K — 12, M = 1 and a — 0.8. The red circles represent the radii of the event and Cauchy horizons, 
while the yellow circles denote the minimal and maximal radius of the orbit, respectively. 




FIG. 11: The three domains Nl to N3 in the Lz-E^-plane for which solutions to the geodesic equation for massless particles 
exist are shown for AI — 1, a = 0.8 and K = 12. The red and blues lines come from the restriction 0{9) > and R{r) > 0, 
respectively. 
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r-e plane 
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FIG. 12: The change of the bound orbit (domain N3) of a massless test particle due to the change of the deficit parameter /?. 
Here Lz = 3.0, E = y/0.90, K = 12, M — 1 and a = 0.8. The red circles represent the radii of the event and Cauchy horizons, 
while the yellow circles denote the minimal and max;imal radius of the orbit, respectively. 




FIG. 13: The change of the cross-over flyby orbit (domain Nl) of a massless test particle due to the change of the deficit 
parameter /3. Here = —0.50, E = \/20, K — 12, M — 1 and a = 0.8. The red circles represent the radii of the event and 
Cauchy horizons, while the yellow circle denotes the minimal radius of the orbit. 
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FIG. 14: The change of the flyby orbit (domain N3) of a massless test particle due to the change of the deficit parameter /3. 
Here L^, = 3, E = \/0.9, K = 12, M = I and a — 0.8. The red circles represent the radii of the event and Cauchy horizons, 
while the yellow circle denotes the minimal radius of the orbit. 



